Kondo effect and RKKY interaction in magnetic trimers 
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Clusters of magnetic atoms lying on a metallic substrate can exhibit Kondo effect due to the ex- 
change interaction with mobile electrons in the metal, which can be observed in STM measurements. 
The same interaction results in the induced RKKY interaction between the magnetic moments in 
the cluster. For different geometries of the cluster, it can enhance or suppress Kondo-correlations. 
We calculate the tunneling density of states for trimer clusters of different geometries and reveal 
the signatures of competition between the Kondo effect and RKKY interaction that can be seen in 
STM experiments. 



Competition of Kondo correlations and RKKY interac- 
tions has a profound effect on physics of magnetic alloys 
and artificially created structures such as quantum dots 
or clusters of magnetic atoms [H, 0, S S] • 

Despite their common origin in the exchange interac- 
tion between localized magnetic impurities and mobile 
electrons, the Kondo correlations and RKKY interac- 
tion lead to quite opposite behavior at low energy scales. 
Kondo correlations lead to the screening of localized mag- 
netic moments by the spins of mobile electrons, resulting 
in the formation of a so-called Kondo cloud around each 
magnetic impurity. The Kondo screening is accompanied 
by a drastic rearrangement in the energy spectrum of mo- 
bile electrons: the Kondo peak in the one particle density 
of states (DoS) appears at the Fermi energy. This in turn 
leads to enhanced scattering of electrons by the screened 
magnetic impurity, which is registered experimentally as 
increase of the resistance of a metal with lowering tem- 
perature, or as increase of conductance through the ar- 
tificial quantum dot P, [1]. In contrast RKKY interac- 
tion leads to magnetic correlations between distant mag- 
netic impurities, thus prohibiting their screening. While 
changing its sign with the distance, RKKY interaction 
can favor ferro- or antiferromagnetic ordering, or result 
in the spin glass order in alloys with random positions 
of magnetic impurities [sl . In the case of dominating an- 
tiferromagnetic RKKY interaction, the Kondo effect is 
suppressed. 

Theoretical investigations of the competition between 
the Kondo effect and RKKY interaction mostly focused 
on the model of two magnetic impurities in a metal 
Inspired by experiments on double quantum dots 
in Kondo regime [4] This model was further extended 
in application to the transport through double quantum 
dot |8"| . There is a general result of theoretical considera- 
tions that strong RKKY interaction binds the two mag- 
netic moments into a singlet state and thus suppresses 
the Kondo effect. 

The interplay between Kondo effect and RKKY inter- 
action becomes even more involved in frustrated systems, 
such as clusters of three magnetic atoms with antifer- 
romagnetic RKKY interaction. Artificial fabrication of 



such clusters deposited on a metallic substrate has be- 
come possible using the scanning tunneling microscopy 
(STM) Hi. Moreover, STM technique allowed a control- 
lable change of the positions of atoms and hence con- 
trolled change of the geometry of the trimer. The ex- 
periments showed clear dependence of low-energy den- 
sity of states in the trimer on its geometry, which can be 
attributed to the competition of the Kondo effect and 
RKKY interactions. Subsequent theoretical investiga- 
tions by means of continuous time Monte Carlo methods 
also showed the change of the density of states with the 
geometry of a triangular atomic cluster, and in partic- 
ular the suppression of the Kondo peak in DoS for the 
antiferromagnetic sign of RKKY coupling 9] . 

Because of the common origin of the Kondo correla- 
tions and RKKY interaction it is important to treat both 
factors on equal footing. In this Letter we present such a 
treatment based on the mean field theory approach sim- 
ilar to the one used for the Kondo effect in the two level 
quantum dot We present results on the local den- 

sity of states (DoS) at positions of magnetic impurities, 
concentrating on the antiferromagnetic (AFM) RKKY 
interaction. There are two qualitatively different regimes 
that are discriminated by the relative strengths of Kondo 
effect and RKKY coupling. In the Kondo regime, which 
is realized by weak RKKY interaction, there are pro- 
nounced Kondo peaks in DoS at positions of the mag- 
netic atoms. The width and the hight of the peaks is 
modified by RKKY interactions and reflects the geomet- 
ric form of the cluster. For an isosceles triangle we find 
a general tendency of the separation of the cluster into 
a single atom and the dimer with growing RKKY inter- 
action. The behavior at strong AFM RKKY coupling 
depends very much on the geometry of the cluster. In 
a strongly asymmetric cluster, two atoms form a singlet 
state with strongly suppressed Kondo effect, while the 
distant atom is singled out and exhibits a Kondo peak 
in the local DoS. In contrast, an almost equilateral tri- 
angle shows a complex evolution of DoS with increasing 
RKKY interaction. At two atoms that are closer to each 
other, the Kondo peak first diminishes and splits, and 
then grows with RKKY coupling. At the same time there 
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is a growth in the local DoS at the third atom. Moreover 
the positions of the peaks in DoS deviates from the Fermi 
level. Our findings are summarized in Figs. [U O 

We consider a cluster of three magnetic atoms each 
having a spin 1/2 placed on a metallic substrate. The 
microscopic description of the system is given by the 
Hamiltonian of s-d model generalized to the case of three 
magnetic impurities 

kcr ja 

+ E ^^^T^ji + E E (^Jk^aak. + h.c. ) (1) 

j j ko- 

Here operators Cja , c]^. relate to the localized d- level in 

the atom j, the operators ako-iakg. relate to the mobile 
fermions in the metallic substrate. The hybridization be- 
tween the localized levels and the extended states in the 
substrate is given by the last term in ([1]). The RKKY 
interaction results from the finite overlap of the extended 
state k with the localized states of different atoms, that 
is 7^ and T^/k ^ for j ^ j'. To take into account 
RKKY interaction on equal footing with the Kondo ef- 
fect, we formally extend the basis of states for mobile 
fermions introducing a three-component field for each ex- 
tended state I kcr) with components 
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interaction term with the help of a Hubbard-Stratonovich 
transformation introducing the hybridization fields as fol- 
lows 



I I' kk' J 



dbj exp 



J + 36, E(*Jk-<f,+$, -*,k' 



kk' 



(6) 



The Schrieffer- Wolff transformation reduced the Hilbert 
space at each magnetic impurity to the two spinful states. 
This allowed us to use a semifermionic representation of 
spin degrees of freedom at each impurity in the path in- 
tegral formulation of partition function l^, In the 
semifermionic representation, the propagator of the cor- 
responding Grassman fields is given by 



i7rT/2 



(7) 



where T denotes the temperature. We made a static 
mean- field ansatz for the decoupling fields, 6J — bj (t 
denotes the imaginary time). Then, after integrating out 
the fermionic variables, the partition function acquires 
the form Z — e^^^^ with the free energy F given by 



where ijk = Tjk/To with To{k) = VETl^jkp- Note that 
this formal extension does not introduce any additional 
degrees of freedom. The free Green's function of the fields 
ipjka in Matsubara frequency space becomes a matrix 
with elements 



[G(iw„,k) 
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The coupling of the fields ^pj and ipji takes into account 
the elementary processes responsible for the RKKY inter- 
action. After performing a Schrieffer- Wolff transforma- 
tion, we obtain the effective Kondo Hamiltonian, which 
can be represented in the form 



= E ^kaL^k. + ^ E E [^Ik'^'J'jk' 



k.k' j 



(4) 

where *jk = [ipji^'i ,ip]i^iV , = {cj],CjiY, and the 
Kondo coupling is given by the standard expression 
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(5) 



In deriving ([5]) we assumed the average tunneling To(k) 
to be independent of the wave vector. To develop a mean 
field treatment of the problem, we write down the parti- 
tion function for the Hamiltonian ([3]), and decouple the 



F = Fo + -E6j -3rETrln] le + DBEGkBi . 

i=i y. k J 

(8) 

Here Fq is an additive constant independent of decou- 
pling fields bj, and the matrices are given by the fol- 
lowing expressions: D(w„) = {iuJn + iirT /2)~^'\.q, B = 
diag(6j) ® 12- The elements of the matrix ^[^Gk(«a;„) 
read 



lUJr, 



Ck 



5aa' = 9jj'ii^n)- 



(9) 

The matrix elements gjji depend on the geometry of the 
magnetic cluster. Taking the positions of atoms as rj, we 
can write i^k = e"*"^""^ /VS. The explicit expression for 
gjj> (itOn) now reads 



3(27r) 
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(10) 

Here rjj' = jr, — rj/|, and we took into account that the 
energy ek depends only on the absolute value k. The ac- 
tual value of gjji depends on the coordinates of magnetic 
atoms in an oscillatory manner and carries the informa- 
tion about RKKY equation. Approximate evaluation of 
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(IT0)1 gives 
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67r2|rj -ryP | 2 




1 - 



+ 1 arctan 
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Here m* denotes the effective mass of an electron, and 
cf denotes the Fermi energy. Perturbative evaluation of 
RKKY interaction in terms of parameters gjj' results in 
the expression 



J2 



(12) 



One can see from (jlip that gjji takes complex values, 
and the relation between its real and imaginary parts 
depends on the distance between the atoms \Tj — Yji\ 
and on the value of u;„. The dependence on the in- 
teratomic distance leads to the changing sign of the 
RKKY interaction. In particular, it follows from (fT2|) 
that real- valued cjjji corresponds to the antiferromagnetic 
RKKY coupling. It is convenient to parameterize gjji as 
gij{iujn) = —isgii(ujn)gjj'{uJn)- In what follows we con- 
sider a cluster in the form of an isosceles triangle with 
AFM RKKY coupling by choosing the explicit form of 
the matrix gjji as follows 



50 m m 
igi go igi 
ig2 igi go 



(13) 



Minimizing the free energy ([8]), we obtain mean field 
equations for the fields bj that for low temperatures can 
be written as 



(14) 

Here Aq is a high-energy cutoff of the order of the band- 
width, which is typical for the Kondo problem, and Xi = 
-|- iA" are the eigenvalues of the matrix BgB. The 
mean field solutions can further be used to determine the 
local one particle density of states in presence of magnetic 
impurities. The single particle Green's function in the 
mean field approximation can be written as 

G(e;k,k') = G°(e;k)(5kk' + G0(e;k)T(e)G°(e;k'). (15) 

Here the matrix elements of the transfer matrix T(e) can 
be calculated as 



(16) 
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FIG. 1: Color online. Induced density of states for a strongly 
asymmetric isosceles triangle. The asymmetry of the triangle 
is fixed by the ration 32/51 = 40,50,70 in the order of the 
diminishing height of the peak. The parameter gi = O.Olgo 
remained fixed. 



where U is the matrix diagonalizing BgB. The density of 
states at a point r is given as p(e,r) = — ilmG''^(e, r, r). 
The expression for the additional density of states in- 
duced by magnetic impurities acquires especially simple 
form at the points where the magnetic atoms are situ- 
ated. It reads 



Aft(e) 



(u-i)„,;a;u 



a;')2 



(17) 



The induced density of states has been calculated for a 
magnetic clusters of different form by solving numeri- 
cally the mean field equations (fT4)) . Taking into account 
the symmetry of the isosceles triangle, we assumed equal 
mean field values of the field bi and 63. 

In the case of a strongly asymmetric isosceles triangle 
there is a clear tendency to formation of magnetic dimer 
with the increase of the asymmetry of the triangle. The 
results for the induced DoS at one of the atoms forming 
a dimer are shown in Fig. [T] The peak in the density 
of states is suppressed with growing antiferromagnetic 
RKKY interaction. A much more complicated behavior 
is found by the cluster in form of an almost equilateral 
triangle. The antiferromagnetic RKKY interaction in- 
troduces frustration of the magnetic state of the cluster, 
which leads to a degenerate ground state. This degener- 
acy is very likely to be broken by a small asymmetry of 
the triangle. Such an asymmetry has been introduced for 
calculations shown in Fig. [21 The calculations are per- 
formed for an isosceles triangle with the form very close 
to the equilateral one. We found that even weak sym- 
metry breaking in presence of frustration leads to large 
changes in local DoS at different atoms. The trimer again 
demonstrates a tendency to split into a dimer and a single 
atom. The density of states however bears signatures of 
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FIG. 2: Color online. Induced density of states for almost 
equilateral triangle (Gi — l.Olgi) for gi/go = 0.5,0.6,0.7. 
The panel a) shows DoS on one of the atome forming a dimer. 
The panel b) shows DoS on a singled out atom. 

strong AFM interaction between the atoms. One can dis- 
tinguish two regimes, the Kondo regime at weak RKKY 
interaction and the regime of frustrated antifcrromag- 
netic cluster at strong interatomic RKKY exchange. In 
the Kondo regime, there is a Kondo peak at the position 
of each atom in the trimer. The Kondo peak is centered 
at the Fermi enegy (e = 0). We found that weak cor- 
relations between magnetic moments that are mediated 
by RKKY interactions lead to the enhancement of the 
Kondo peak (not shown). 

At strong RKKY interaction the regime of antiferro- 
magnetic cluster is realized. DoS reflects the separation 
of the equilateral triangle into a dimer and a singled out 
atom. This is indicated by a strong drop in DoS on dimer- 
ized atoms accompanied by a splitting of the Kondo peak 
(see the curve for gi/go = 0.6 in the panel a) while there 
is a pronounced Kondo peak at the position of the sep- 
arated atom (panel b). The splitting of the Kondo peak 
has also been observed in numerical simulations [!| . With 
further growth of RKKY interaction, the splitting of the 
Kondo peak at a dimer atom becomes asymmetric in en- 
ergy, the peak in DoS shifts from the position of the Fermi 
level (the curve gi/ go = 0.6 in the panel a). DoS on the 
singled out atom is shown in the panel b) to Fig. [21 The 
Kondo peak at the singled out atom becomes broader 
and lower with growing RKKY coupling. We attribute 
this peculiar behavior to the strong frustration in the 
equilateral triangular cluster with AFM nearest neigh- 
bor interaction. The shift of the Kondo peak from the 
position of the Fermi level (e = 0) may indicate an ap- 
pearance of a new phase for the triangular cluster which 
might demonstrate a non- Fermi-liquid behavior [14]. 

In conclusion, we proposed a theoretical treatment that 
takes into account Kondo correlations and RKKY in- 



teraction between localized magnetic moments on equal 
footing. We applied this approach to triangular mag- 
netic clusters and investigated the dependence of the lo- 
cal density of states on the geometry of the cluster. We 
found two qualitatively different regimes. In the regime 
of weak antiferromagnetic RKKY interaction, the RKKY 
coupling enhances Kondo peak in DoS at positions of 
magnetic atoms. In the opposite regime of strong RKKY 
interaction, the Kondo effect is strongly suppressed, and 
DoS is essentially featureless. 

We believe that our approach can further be applied 
to the magnetic clusters in the mixed valence regime and 
to the clusters with larger number of atoms. The form of 
the mean field equations and of the local density of 
states ([T7| remains unchanged, while the matrix g (fT3)) 
for a cluster of N atoms becomes A^-dimensional. Di- 
agonalizing the matrix BgB thus becomes more difficult 
but should be numerically accomplished without severe 
problems. 
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